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Tenth-Order QED Contribution to the Electron g — 2 
and an Improved Value of the Fine Structure Constant 
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This paper presents the complete QED contribution to the electron g~2 up to the tenth order. 
With the help of the automatic code generator, we have evaluated all 12672 diagrams of the tenth- 
order diagrams and obtained 9.16 (58)(Q/7r)^. We have also improved the eighth-order contribution 
obtaining —1.9097 (20)(a/7r)^, which includes the mass-dependent contributions. These results lead 
to ae(theory) = 1 159 652 181.78 (77) x 10^^^. The improved value of the fine-structure constant 
= 137.035 999 174 (35) [0.25ppb] is also derived from the theory and measurement of Oe. 

PACS numbers: 13.40.Em, 14.60.Cd, 06.20.Jr, 12.20.Ds 



The anomalous magnetic moment = {g — 2)/2 of 
the electron has played the central role in testing the 
validity of quantum electrodynamics (QED) as well as 
the standard model of the elementary particles. On the 
experimental side the measurement of ae by the Harvard 
group has reached the astonishing precision P, Q : 

ae(HV) = 1 159 652 180.73 (0.28) x lO'^^ [0.24ppb](l) 

In the standard model the contribution to Og 
comes from three types of interactions, electromagnetic, 
hadronic, and electroweak: 

fle = ae(QED) + fle (hadronic) -I- ag (electroweak) . (2) 

The QED contribution can be evaluated by the pertur- 
bative expansion in a/n: 



integration [1, They are summarized as: 



a(^) = 0.5, 

A^^^ = -0.328 478 965 579 193 



1.181 241 456 .. 



= -1.9106 (20), 



A 



(10) 



9.16 (58) 



(5) 
(6) 



(8) 

The A\ is obtained from 891 Feynman diagrams clas- 
sified into 13 gauge-invariant subsets (see Fig. [T|). The 
value = -1.9144 (35) in Q was confirmed by the 
new calculation and replaced by the updated value ([S|). 

The A^^^-* receives the contribution from 12672 di- 
agrams classified into 32 gauge-invariant subsets (see 
Fig. [2]). The results of 31 gauge-invariant subsets have 
been published [IMi- The remaining set. Set V, con- 
sists of 6354 diagrams, which are more than half of all 
tenth-order diagrams. However, we have managed to 
evaluate it [20| with a precision which leads to theory 
more accurate than that of the measurement ([l} : 

A^^°^[Set V] = 10.092 (570). 



ae(QED) = £(^)"a(2n)^ 



(3) 



ra=l 



(7) 

where a?"^ is finite due to the renormalizability of QED ^^ding data of all 32 gauge-invariant subsets, we are now 



and may be written in general as 



-I- A^^"''(me/m^,me/mr) (4) 



to show the mass-dependence explicitly. We use the lat- 
est values of the electron- muon mass ratio me/m^ = 
4.836 331 66 (12) x 10~^ and the electron-tau mass ratio 
me/mr = 2.875 92 (26) x 10"^ 1]. 

The first three terms of ^^^"^ are known analytically 
[3-0|! while A^^'^ and A^^^^ are known only by numerical 



able to obtain the complete value of A^^"' as in ([6|) , which 



replaces the crude estimate A^^"^ = 0.0(4.6) SmUll. 

The mass-dependent terms Ai and A^ of the fourth 
and sixth orders are known [231428 1 and re-evaluated us- 
ing the updated mass ratios [3|, 



A':^\mjm^,) ^ 5.197 386 67 (26) x 10" 
A':^\melmr) = 1.837 98 (34) x 10"^ 



Af\me/mi^, 



Af\me/mr) 
Af'' {me/m^,me/mr 



7.373 941 55 (27) x 10-^ 
6.583 (11) X 10"^ 

0.1909 (1) X 10"^^ 



(8) 



Except for A^^ all are known analytically so that the 
uncertainties come only from fermion-mass ratios. 

The mass-dependent terms of the eighth order and the 
muon con tributio n to the tenth order are numerically 
evaluated 



11-191. Our new results are summarized as 



A^^^\me/m^) = 9.222 (66) x lO"'^, 
A^^\me/mr) ^ 8.24 (12) x 10"^, 
A'f\me/mf„me/mr) = 7.465 (18) x 10~^ 
A^^^°\me/m^) = -0.003 82 (39). 



(9) 



The hadronic contribution to ae is summarized in 
Ref. 0]. The leading order [2^ and next-to- leading or- 
der (NLO) [i^l contributions of the hadronic vacuum- 
polarization (v. p.) as well as the hadronic light- by- light- 
scattering {l-l) term [3l| are given as 



ae(had. v.p.) = 1.875 (18) x 10"^^ 
ae(NLO had. v.p.) = -0.225 (5) x 10"^^ 
ae(had. l-l) = 0.035 (10) x 10"^^ 



(10) 



At present no direct evaluation of the two-loop elec- 
troweak effect is available. The best estimate is the one 
obtained by scaling down from the electroweak effect on 



ae(weak) = 0.0297 (5) x 10" 



(11) 



To compare the theoretical prediction with the mea- 
surement ([T]), we need the value of the fine-structure con- 
stant a determined by a method independent oi g — 2 . 
The best a available at present is the one obtained from 
the measurement of h/m^^ 36|, combined with the very 
precisely known Rydberg constant and mjib/me Q : 

a"i(Rbl0) = 137.035 999 049 (90) [0.66ppb]. (12) 

With this a the theoretical prediction of becomes 



ae(theory) =1 159 652 181.78 (6)(4)(3)(77) 



X 10" 



[0.67ppb], 



(13) 



where the first, second, third, and fourth uncertainties 
come from the eighth-order term ([S]), the tenth-order 
term ([6|), the hadronic corrections (ITO|) . and the fine- 
structure constant p^ . respectively. This is in good 
agreement with the experiment ([Ij: 

ae(HV) - ae(theory) = -1.06 (0.82) x 10"^^ (14) 

More rigorous comparision between experiment and the- 
ory is hindered by the uncertainty of a~^(Rb) in ([T2|). 
Note that the sum 1.685(21) x 10"^^ of the hadronic 
contributions (|10|) is now larger than Eq. (fT4|). It is 
thus desirable to reexamine and update the values of the 
hadronic contributions. 



\ 
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FIG. 1. Typical vertex diagrams representing 13 gauge- 
invariant subsets contributing to the eighth-order lepton <?— 2. 



The equation (|T3|) shows clearly that the largest source 
of uncertainty is the fine-structure constant (1121) . To put 
it differently, it means that a non-QED a, even the best 
one available at present, is too crude to test QED to 
the extent achieved by the theory and measurement of 
flg. Thus it makes more sense to test QED by an alter- 
native approach, namely, compare a^^(RblO) with a^^ 
obtained from theory and measurement of Ue- This leads 
to 

a-i(ae) = 137.035 999 1736 (68)(46)(26)(331) [0.25ppb], 

. 

where the first, second, third, and fourth uncertainties 
come from the eighth-order and the tenth-order QED 
terms, the hadronic and electroweak terms, and the mea- 
surement of Oe (HV) in (m , respectively. The uncertainty 
due to theory has been improved by a factor 4.5 com- 
pared with the previous one ^22;]. 

Let us now discuss the eighth- and tenth-order calcu- 
lations in more details. The 13 gauge-invariant groups of 
the eighth order were numerically evaluated by VEGAS 
[39I and published 0, [l3| • As an independent check, we 
built all programs of the 12 groups from scratch with 
the help of automatic code generator GENCODeA, except 
for Group IV(d) which had already been calculated by 
two different methods The new values of the mass- 
independent contributions of all 12 groups are consistent 
with the old values. We have thus statistically combined 
two values and listed the results in Table I. Since the va- 
lidity of the new programs were confirmed in this way, 
we used the new programs as well as the old programs to 
evaluate the mass-dependent terms A'^^ and ^3^''. 

Group V deserves a particular attention which consists 
of 518 vertex diagrams and is the source of the largest un- 
certainty of ai*'' . The programs generated by GENCODeA^ 
have been evaluated with intense numerical work which 
led to —2.173 77 (235). This is consistent with the value 
in H, -2.179 16 (343). The combined value is 



[Group V] = -2.175 50 (194) 



(16) 



This improvement results in about 40% reduction of the 
uncertainty of the eighth-order term. 

The tenth-order contribution comes from 32 gauge- 
invariant subsets (see Fig. [2]). The FORTRAN programs 
of integrals of 15 subsets I(a-f), II(a,b), 11(f), Vl(a-c), 
VI(e,f), and VI(i) are straightforward and obtained by a 



3 



TABLE I. The eighth-order QED contribution from 13 gauge-invariant groups to electorn (7—2. The values with a superscript 
a, b, or c are quoted from Refs.js^], @|, or [s^, respectively. Uf shows the number of vertex diagrams contributing to A^^\ 
Other values are obtained from evaluation of new programs. The mass-dependence of Ag*^ is A^^\me/m^,me/mT). 



group 


Uf 




A^^\mjm^) X 10^ 


Ai^^^mjuir) X 10^ 


4«) X 10^ 


1(a) 


1 


0.000 876 865 ■ ■ 


0.000 226 456 (14) 


0.000 080 233 (5) 


0.000 


Oil 994 


1(b) 


6 


0.015 325 20 (37) 


0.001 704 139 (76) 


0.000 602 805 (26) 


0.000 


014 097 


1(c) 


3 


0.011 130 8 (9)'' 


0.011 007 2 (15) 


0.006 981 9 (12) 


0.172 


860 (9^ \ 


1(d) 


15 


0.049 514 8 (38) 


0.002 472 5 (7) 


0.087 44 (1) 







11(a) 


36 


-0.420 476 (11) 


-0.086 446 (9) 


-0.045 648 (7) 







11(b) 


6 


-0.027 674 89 (74) 


-0.039 000 3 (27) 


-0.030 393 7 (42) 


-0.458 


968 (17) 


11(c) 


12 


-0.073 445 8 (54) 


-0.095 097 (24) 


-0.071 697 (25) 


-1.189 


69 (67) 


III 


150 


1.417 637 (67) 


0.817 92 (95) 


0.6061 (12) 







IV(a) 


18 


0.598 838 (19) 


0.635 83 (44) 


0.451 17 (69) 


8.941 


(17) 


IV(b) 


60 


0.822 36 (13) 


0.041 05 (93) 


0.014 31 (95) 







IV(c) 


48 


-1.138 52 (20) 


-0.1897 (64) 


-0.102 (11) 







IV(d) 


18 


-0.990 72 (10)"= 


-0.1778 (12) 


-0.0927 (13) 







V 


518 


-2.1755 (20) 
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FIG. 2. Typical self-energy-like diagrams representing 32 
gauge- invariant subsets contributing to the tenth-order lepton 
g — 2. Solid lines represent lepton lines propagating in a weak 
magnetic field. 

slight modification of programs for the eighth-order dia- 
grams. Together with the results of subsets VI(j,k), the 
contributions from 17 subsets to A^^*^-* were evaluated 
and published We recalculated all 17 subsets once 
more from scratch and found that the results of 1(d), 1(f), 
11(a), 11(b), and VI(c) in [l3| were incorrect. Although 
the constructed integrals for the first four subsets are free 
from errors, they did not include the finite renormaliza- 
tion terms in the last step of the calculation. The value 
of the subset VI(c) was a typo. The corrected values are 
listed in Table II. 

Other subsets are far more difficult to handle. Thus 
we developed and utilized the code-generating algorithm 
GENCODeA^ which carries out all steps automatically. 



including subtraction of ultraviolet and infrared diver- 
gences |40|. By GENCODeA^ and its modifications for 
handling vacuum-polarization loops and light-by-light- 
scattering loops, we have obtained FORTRAN programs 
for 12 more subsets [ll[li-[li|. The subsets III(c) and 
I(j), which involve one(two) light- by- light scattering sub- 
diagram(s) internally, were handled manually |lll [l9| . 
The subset 11(e), which contain a sixth-order light-by- 
light-scattering subdiagarm internally, was handled by 
an automation procedure [l^. At least two independent 
codes for non-automated programs were written by dif- 
ferent members of our collaboration in order to minimize 
human errors. 

All integrals were numerically evaluated by VEGAS 
(sof . For some diagrams of the sets IV and V that contain 
cancellation of linear IR divergence within a diagram, we 
used the quadruple-precision arithmetics to avoid possi- 
ble round-off errors of numerical calculations. The con- 
tribution of the tau-particle loop to is negligible at 
present. Thus the sum of (O and (jH]) gives effectively the 
total tenth-order QED contribution to Og. 

This work is supported in part by the JSPS 
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(C)23540331. T. K.'s work is supported in part by the 
U. S. National Science Foundation under Grant NSF- 
PHY-0757868. T. K. thanks RIKEN for the hospitality 
extended to him while a part of this work was carried 
out. Numerical calculations are conducted on RSCC and 
Rice supercomputer systems at RIKEN. 
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